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ABSTRACT 

An accurate geometric distortion solution for the Hubble Space Telescope 
UVIS-channel of Wide Field Camera 3 is the first step towards its use for high 
precision astrometry. In this work we present an average correction that enables 
a relative astrometric accuracy of ~1 mas (in each axis for well exposed stars) 
in three broad-band ultraviolet filters (F225W, F275W, and F336W). More data 
and a better understanding of the instrument are required to constrain the solu- 
tion to a higher level of accuracy. 



1 Visiting PhD Student at STScI under the "2008 graduate research assistantship" program. 
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1. Introduction, Data set, Measurements 

The accuracy and the stability of the geometric distortion^ (GD) correction of an in- 
strument is at the basis of its use for high precision astrometry. The particularly advanta- 
geous conditions of the Hubble Space Telescope (HST) observatory make it ideal for imaging- 
astrometry of (faint) point sources. The point-spread functions (PSFs) are not only sharp and 
(essentially) close to the diffraction limit -which directly results in high precision positioning- 
but also the observations are not plagued by atmospheric effects (such as differential refrac- 
tion, image motion, differential chromatic refraction, etc), which severely limit ground-based 
astrometry. In addition to this, HST observations do not suffer from gravity- induced flexures 
on the structures of the telescope (and camera), which add (relatively) large instabilities in 
the GD of ground-based images, and make its corrections more uncertain. 

Last May 14, the brand-new Wide Field Camera 3 (WFC3) was successfully installed 
during the Hubble Servicing Mission 4 (SM4, May 12-24 2009). After a period of intense 
testing, fine-tuning, and basic calibration, last September 9th, 2009, the first calibration- 
and science-demonstration images were finally made public. 

Our group is active in bringing HST to the state of the art of its astrometric capabilities, 
that we used for a number of scientific applications (e.g. from King et al. 1998, to Bedin et al. 
2009, first and last accepted papers). Now that the "old" ACS/WFC is successfully repaired, 
and that the new instruments are installed, our first step is to extend our astrometric tools to 
the new instruments (and to monitor the old ones). This paper is focused on the geometric 
distortion correction of the UV / Optical (UVIS) channel of the WFC3. Since the results of 
these efforts might have some immediate public utility (e.g. relative astrometry in general, 
stacking of images, UV-identification of X-counterparts such as pulsars and CVs in globular 
clusters, etc.), we made our results available to the WFC3/UVIS user- community. 

We immediately focused our attention on a deep UV-survey of the core of the Galactic 
globular cluster uj Centauri (NGC 5139), where some well dithered images were collected. 
The dense -and relatively flat- stellar field makes the calibration particularly suitable for 
deriving and monitoring the GD on a relatively small spatial scale. In addition, while most 
of the efforts to derive a GD correction will be concentrated on relatively redder filters, we 
undertook a study to determine the GD solutions of the three bluest broad-band filters (with 
the exception of F218W): F225W, F275, and F336W. 



1 A specification is needed. With the term "geometric distortion", or GD, which we will use hereafter, 
we are lumping together several effects under the same term: the optical field-angle distortion introduced 
by camera optics, light-path deviations caused by the filters, non-flat CCDs, alignment errors of CCDs on 
the focal plane, etc.. 
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The WFC3/UVIS layout is almost indistinguishable from that of ACS/WFCB two E2V 
thinned, backside illuminated and UV optimized 2kx4k CCDs contiguous on the long side 
of the chip, and covering a field of view (FoV) of ~ 160x160 arcsec 2 . The u Cen data set 
used here consists of 9x350 s exposures in each of the filters F225W, F275W, and F336W. 
The nine images follow a squared 3x3 dither-pattern with a step of about 40 arcsec (i.e. 
~1000 pixels), and were all collected on July 15, 2009. 

We downloaded the standard pipe-line reduced FLT files from the archive. The FLT 
images are de-biased and flat-field corrected, but no pixel-resampling is performed on them. 
The FLT files are multi-extension fits (MEF) on which the first slot contains the image of 
what -hereafter- we will call chip 1 (or simply [1]). The second chip, instead, is stored in 
the fourth slot of the MEF, and we will refer to it as chip 2 (or [2]). [Note that others might 
choose a different notation] . Our GD corrections refer to the raw pixel coordinates of these 
images. 

The fluxes and positions were obtained from a code mostly based on the software 
img2xym_WFI by Anderson et al. (2006). This is essentially a spatially variable PSF-fitting 
method. We were pleased to see that for the WFC3/UVIS images of this data set the 
PSFs were only marginally undersampled. Left panel of Figure [1] shows a preliminary color- 
magnitude diagram in the three filters for the bright stars in the WFC3/UVIS data set. In 
a future paper of this series we will discuss the PSF, its spatial variation and stability, as 
well as L-flat9§, pixel-area corrections, and recipes for deep photometry in stacked- images. 



2 http://www.stsci.edu/hst/wfc3 

3 Residual low-frequency flat-field structure (L-flat) cannot be accurately determined from ground-based 
calibration data or internal lamp exposures. L-flats need to be determined from on-orbit science data, for 
example from multiple observations of stellar fields with different pointings and roll angles (van der Marel 
2003). 
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Fig. 1.— Left: Preliminary color- magnitude diagram of the bright stars in the new 
WFC3/UVIS data set (fluxes are neither pixel-area- nor L-flat-corrected). Right: Color- 
magnitude diagram of the stars in our ACS/WFC master frame (from Villanova et al. 2007). 
Both plots are in instrumental magnitudes. 
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Fig. 2. — Left: Color- magnitude diagrams of the high S/N stars in common with the 
ACS/WFC master frame (F435W), actually used to derive the geometric distortion cor- 
rection, for each of the WFC3/UVIS filters (F225W, F275W, and F336W). 
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2. 



The Geometric Distortion Solutions 



The most straightforward way to solve for the GD would be to observe a field where there 
is a priori knowledge of the positions of all the stars in a distortion-free reference frame. [A 
distortion-free reference frame is a system that can be transformed into any other distortion- 



itself immediately as the residuals between the observed relative positions of stars and the 
ones predicted by the distortion- free frame (on the basis of a conformal transformation). 

Thankfully, we possess such an "astrometric flat-field", moreover with the right mag- 
nitude interval, source density, and accuracy. This reference frame is the mosaic of 3x3 
Wide Field Channel (WFC) of the Advanced Camera for Surveys (ACS) fields collected - 
at the end of June 2002- under the program GO-9442 (PI: Cool) reduced by Jay Anderson 
and published in Villanova et al. (2007). This reference frame was obtained from a total of 
18 short and 90 long ACS/WFC exposures in the filters: F435W, F625W, and F658N (see 
Villanova et al. 2007 for details). The entire field covers an area of 10x10 arcmin 2 , and can 
be considered distortion-free at the ~0.5 mas level. The catalog contains more than 2 million 
sources, and we will refer to it as master frame, and to the coordinates of each i-source in it 
with the notation (X J master ,y i master ). A color-magnitude diagram for the stars in the master 
frame is shown on the right panel of Figure [TJ 

To derive the WFC3/UVIS GD corrections we closely follow the procedures described 
in Anderson & King (2003, hereafter AK03) used to correct the GD for each of the four 
detectors of WFPC2. We represent our solution with a third-order polynomial, which is able 
to provide our final GD correction to the ~0.025 pixel level in each coordinate (~1 mas). 
Higher orders proved to be unnecessary. 

Having a separate solution for each chip, rather than one that uses a common center 
of the distortion in the FoV, allows a better handle of potential individual detector effects 
(such as a different relative tilt of the chip surfaces, relative motions, etc.). We adopted as 
the center of our solution, for each chip, the point (x ,y )k=i,2 = (2048,1025) [in the raw 
pixel coordinates, to which we will refer hereafter as (x, y)\. 

For each i-star of the master list, in each /c-chip of each j-MEF-file, the distortion 
corrected position (x corr , y corr ) is the observed position (x,y) plus the distortion correction 



free frame by means 



of conformal transformation^].] Geometric distortion would then show 




A conformal transformation between two catalogs of positions is a four-parameter linear transformation, 
specifically: rigid shifts in the two coordinates, one rotation, and one change of scale, i.e. the shape is 
preserved. 
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(5x,5y): 



f X i,j,k X i,j,k + $ x i,j,k( x i,j,k,yi,j,k) 
I Vi°j,k = Vi,j,k + $yi,j,k(Xi,j,k, 2/i,j,fc), 

where and y^-^ are the normalized positions, defined as: 




Normalized positions make it easier to recognize the magnitude of the contribution given by 
each solution term, and their numerical round-off. 

The final GD correction for each star, in each chip/image, is given by these two third- 
order polynomials (we omitted here i, j, k indexes for simplicity): 



Our GD solution is thus fully characterized by 18 coefficients: ai,...,ag, bi, . . . , bg. 
However, as done in AK03, we constrained the solution so that, at the center of the chip, 
it will have its rr-scale equal to the one at the location (x ,y ), and the corrected axis y COTT 
has to be aligned with its y-axis at the location (x Q ,y Q ). This is obtained by imposing 
ai t k = and a 2) k = 0. Since the detector axes do not necessarily have the same scales nor are 
perpendicular to each other, b\ t u and &2,fc must be free to assume whatever values fit best. 
Therefore, we have to compute in fact only 16 coefficients (for each chip) to derive our GD 
solution. 

Each i-star in the master frame is conformally transformed into each /c-chip/j-image, 
and cross-identified with the closest source. We indicate such transformed positions with 
(X™ aster ) Tj '- fc and (^ master ) T j,fc. Each of such cross-identifications, when available (of course 
not all the red sources in the master list were available in the WFC3 UV-filters) , generates 
a pair of positional residuals: 



which reflect the residuals in the GD (with the opposite sign), and depend on where the 
i-star fell on the k/j -chip/image (plus random deviations due to non-perfect PSF-fitting and 
photon noise). Note that our calibration process is an iterative procedure, and that necessar- 
ily, at the first iteration, we have to impose (x COTT ,y COTT )ij jk = {x^y)^^. In each chip/image 



5x = a\x + ci2y + a^x 2 + a^xy + a^y 2 + a^x 3 + ajx 2 y + a$xy 2 + agy 
5y = b\x + b 2 y + b 3 x 2 + b^xy + b 5 y 2 + 6 6 5 3 + b 7 x 2 y + b 8 xy 2 + bgy 3 . 
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we have typically ~5500 high-signal unsaturated stars in common with the master frame, 
leading to a total of ~50 000 residual pairs per chip. [A color- magnitude diagram of the stars 
actually used to compute the GD solution is shown, for each filter, in FigO] 

For each chip, these residuals were then collected into a look-up table made up of 37 x 19 
elements, each related to a region of 110 x 110 pixels. We chose this particular grid setup 
because it offers the best compromise between the need for an adequate number of grid 
points to model the GD, and an adequate sampling of each grid element, containing at least 
60 pairs of residuals. For each grid element, we computed a set of five 3cr-clipped quantities: 
x m ,n,k, y m>n ,k' ^ x m,n,k, Ay mnk , and Pm,n,k\ where x m ,n,k and y m ^ k are the averaged positions 
of all the stars within the grid element (m, n) of the /c-chip, Ax mtn ^ and Ay mnk are the 
average residuals, and P m , n ,k is the number of stars that were used to calculate the previous 
quantities. These P m ,n,k will also serve in associating a weight to the grid cells when we fit 
the polynomial coefficients. 

To obtain the 16 coefficients describing the two polynomials (a qt k with q = 3, . . . ,9, and 
b qyk with q = 1, . . . , 9), which represent our GD solution in each chip, we perform a linear 
least-square fit of the iV = mxn = 37 x 19 = 703 data points. Thus, for each chip, we can 
compute the average distortion correction in each cell ( Sx Ptk , Sy pk ) with N relations of the 



(where t\ iPik = x Pi k, t 2y k = y p ^ ■ ■ ■ > ^9,fc = V 3 p,k)y an d where the 16 unknown quantities - 
a q> k and b q ^ — are our fitting parameters (16 for each chip) . 

In order to solve for a q ^ and h q ^, we formed, for each chip, one 9x9 matrix A4 k and 
two 9x1 column vectors V a ,k and Vb,k'- 



form: 





V a , k 




V b>k 



( J2 P Pp,kti,p,k^y P ,k \ 

^2 P Pp,kh, P ,kAy P}k 



p 1 P b 9,p,k 



V Pp,kts,P,k^ X P,k / 



\ Y,pPp,kt9,p,k^y p ,k / 
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The solution is given by two 9x1 column vectors Ak and Bk, containing the best fitting 
values for a q ^ and b qj k, obtained as: 



\ Ct9,k J 



>%k 



Mk- l V h .k. 



\ h,k J 



With the first set of calculated coefficients a q ^ and b q ^ we computed the corrections 
fixi ,j,k and Syi ; j t k to be applied to each i-star of the /c-chip in each j-M EF file, but actually we 
corrected the positions only by half of the recommended values, to guarantee convergence. 
With the new improved star positions, we start-over and re-calculated new residuals. The 
procedure is iterated until the difference in the average correction from one iteration to the 
following one — for each grid point — became smaller than 0.001 pixels. Convergence was 
reached after ~ 100 iterations. The coefficients of the final GD solutions for the two chips, 
and for the three different filters, are given in Table [H 

In Figure [3] we show for the intermediate filter F275W the total residuals of uncorrected 
star positions vs. the predicted positions of the master frame, which is representative of our 
GD solutions. For each chip, we plot the 37 x 19 cells used to model the GD, each with its 
distortion vector magnified (by a factor of x8 in x, and by a factor of xl.5 in y). Residual 
vectors go from the average position of the stars belonging to each grid cell (x, y) to the 
corrected one (x corr , y covv ). Side panels show the overall trends of the individual residuals 
Sx, 5y along x and y directions (where for clarity we plot only a 40% sub-sample, randomly 
selected). It immediately strikes the large linear terms in y, reaching up to ~140 pixels. 

In Figure H] we show, in the same way, the remaining residuals after our GD solution is 
applied. This time we magnified the distortion vectors by a factor xl500 in both axes. 

At this point it is very interesting to examine the r.m.s. of these remaining residuals, 
that show a rather large ~0.15 pixels dispersion. We will see, in the following, that this 
dispersion can be interpreted as the effect of the internal motions of the cluster stars on the 
time baseline of ~7 years between the ACS/WFC observations of the reference frame, and 
the new WFC3/UVIS data set. Indeed, assuming 1) a distance of 4.7 kpc for u> Cen (van der 
Marel & Anderson 2009), 2) an internal velocity dispersion of ~18 km s _1 in our fields (van 
de Ven et al. 2006), and 3) an isotropic velocity distribution for stars, we would expect to 
observe in ~7 years a dispersion of the displacements of ~5.5 mas. This dispersion, assuming 
a pixel scale of ~40 mas for WFC3/UVIS, corresponds to a displacement of ~0.14 pixels 
(also in good agreement with the recent measurements by Anderson & van der Marel 2009). 



To show this more clearly we intercompare the average positions of the nine WFC3/UVIS 
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corrected catalogs in the filter F275W, with those in the F336W, for the stars in common 
between the two filters. All these images were taken at the same epoch, and so positions 
of stars are not affected by internal motion effects. The 1-dimension dispersion should 
reflect our accuracy, and indeed the observed residuals -in this case- have a dispersion 
of ~0.025 pixels (i.e. ~1 mas). Figure [5] illustrates the two situations. On the left-panel 
we plot the displacements between the ACS/WFC epoch of the reference frame and the 
new WFC3/UVIS epoch, while on the right-panel we show the displacements between our 
corrected position in filter F275W and the corrected positions in F336W. On the left-panel 
the internal motions of to Cen dominate the dispersion, while on the right-panel, there are 
no internal motions at all, and what we are left with are our errors only. 

Unfortunately the WFC3/UVIS images were either not enough, or not well dithered 
to perform a pure auto-calibration, and we had to use the ACS/WFC reference frame. 
Nevertheless, even a dispersion of 0.15 pixel within a given cell should be reduced to less 
than 0.02 pixels if averaged over more than 60 residuals. And this should be regarded as an 
upper limit, since we are using 703 grid points to constrain 16 parameters. 

For this reason, the estimated 0.025 pixel accuracy is larger than we would have ex- 
pected. We can not exclude that these residuals could be due to a deviation from an isotropic 
distribution of the internal motion of u Cen (i.e. at the level of < 3 km s _1 ), or simply by un- 
expectedly large errors in the adopted ACS/WFC astrometric flat-field (the master frame). 
Another possibility is that there could be some unexpected (and so far undetected) man- 
ufacturing artifact in the WFC3/UVIS detectors which could affect the positions (such as 
those identified on WFPC2 CCDs, and characterized by Anderson & King 1999, or those of 
ACS/WFC found by Anderson 2002). Finally, it could simply be a higher frequency spatial 
variation which can not be properly represented by a polynomial of a reasonable order, but 
rather by a residual table as done in Anderson (2006) for ACS/WFC. Surely, more data 
are needed to further improve the GD solutions presented in this work, as well as a better 
time-baseline for the understanding of its variations. We want to end this section by pointing 
out that the detection of the internal motions among the stars of a Galactic globular cluster 
is a rather challenging measurement, and it could well be one of the best demonstrations of 
the goodness of our derived geometric distortion solutions. 
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Table 1: The coefficients of the third-order polynomial for each chip and filter. 



Term(w) Polyn. a WiK b w ^ a w ^ 2 \ K,[2\ 



WFC3/UVIS filter F225W 



1 


X 


0.000 


129.230 


0.000 


140.270 


2 


y 


0.000 


1.935 


0.000 


-4.215 


3 


z 2 


12.120 


0.591 


12.021 


0.773 


4 


xy 


-6.279 


5.553 


-6.057 


5.496 


5 


y 2 


0.064 


-3.227 


0.001 


-3.058 


G 


z 3 


0.176 


0.029 


0.149 


0.156 


1-7 

7 


Z 2 y 


-0.057 


0.033 


0.022 


—0.009 


8 


xy 2 


0.004 


-0.041 


0.061 


-0.026 


9 


y 3 


0.035 


-0.023 


0.032 


0.028 




WFC3/UVIS 


FILTER F275W 




1 


X 


0.000 


129.270 


0.000 


140.285 


2 


y 


0.000 


1.925 


0.000 


-4.221 


3 


Z 2 


12.102 


0.581 


12.016 


0.781 


4 


xy 


-6.284 


5.547 


-6.040 


5.493 


5 


y 2 


0.061 


-3.241 


0.001 


-3.048 


6 


X 3 


0.178 


0.033 


0.144 


0.163 


7 


Z 2 y 


-0.056 


0.054 


0.026 


0.007 


8 


xy 2 


0.005 


-0.041 


0.051 


-0.025 


9 


y 3 


0.033 


-0.012 


0.032 


0.020 




WFC3/UVIS 


FILTER F336W 




1 


X 


0.000 


129.438 


0.000 


140.315 


2 


y 


0.000 


1.786 


0.000 


-4.322 


3 


x 2 


12.091 


0.676 


11.994 


0.672 


4 


xy 


-6.188 


5.565 


-6.135 


5.476 


5 


y 2 


0.065 


-3.155 


0.004 


-3.152 


6 


z 3 


-0.062 


0.004 


-0.151 


0.189 


7 


x 2 y 


-0.097 


0.034 


0.074 


-0.027 


8 


xy 2 


0.016 


-0.061 


0.040 


0.005 


9 


y 3 


0.033 


0.016 


0.033 


0.014 



-12- 



WFC3/UVIS-F275W distortion map 



13 




a 

a 
fa 



[WFC3 pix] 



Fig. 3. — Predicted vs. uncorrected positions. The size of the residual vectors is magnified 
by a factor of x8 in i and xl.5 in y. For each chip we plot also individual residuals as 
function of x and y axes. Units are expressed WFC3/UVIS pixels in the reference positions 
(x , y a ). For clarity only a random 40% of the residuals is plotted. 
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WFC3/UVIS-F275W residual map 
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Fig. 4. — Same as Fig. [3] after the correction was applied. The size of the residuals is now 
magnified by a factor of 1500. 
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ACS/WFC(2002.49) vs. 
WFC3/UVIS(2009.54) 



WFC3/UVIS(F336W) vs. 
WFC3/UVIS(F275W) 
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Fig. 5. — Left: Vector-point diagram of displacements for the stars in common between the 
ACS/WFC epoch of the master catalog, and the average of the corrected WFC3/UVIS new 
data in filter F275W. The internal motions of to Cen dominate the observed dispersions, but 
do not prevent a GD solution accurate to ~0.025 WFC3 pixels. Right: Vector-point diagram 
of displacements for the stars in common between the corrected WFC3/UVIS data in filter 
F336W, and those corrected for F275W. The images are collected at the same epoch, and 
no sizable internal motions are present. In this case, the dispersion reflects our errors. The 
red circles, in both panels, indicate the 1-dimensional dispersion of the residuals, and all 
quantities are expressed in units of WFC3/UVIS pixels. 
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Table 2: Interchip transformation parameters. Chip [1] parameters are indicated only for 
clarity. For chip [2] formal errors are given. 



fc-chip a [k] /a {1] 9 [k] -9 {1] (x [ * ] )<^ rf 1 )^ 
[number] ["] [pixel] [pixel] 



[1] 1.00000 0.0000 2048.00 1025.00 

[2] 1.00595 0.0654 2046.00 3098.34 
±2/100 000 ±0.001 ±0.03 ±0.03 



3. Interchip transformations 



For many applications it would be useful to transform the GD corrected positions of 
each chip into a common distortion-free reference frame. We could then, simply conformally 
transform the corrected positions of chip [k] into the distortion corrected positions of chip [1], 
using the following relations: 



/ ™corr \ 
X [l] 




cos(0[i] - 


- Q[k]) 


sin(0[i] - 9[ k ] 




Ot[k] 


_ -sin(0[i] 


-0 { k)) 


cos(6»[i] - B[k) 



x [k] 

, con 

V[k) 



- 2048 

- 1025 



+ 



(rrl^l Ncorr 
{■lo )^ 

Wl \corr 
3 ' [1] 



(y 



where — following the formalism in AK03 — we indicate the scale factor as as ce[k], 



the 

in the corrected 

reference system of chip [1] as (x\ i )m rr and {y^)[°™ ■ Of course, for k = 1, we end up with 
the identity. The values of the interchip transformation parameters are given in Table [21 and 
shown for individual images in Figure El 



orientation angle with 9^, and the positions of the center of the chip (x , y 
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1.0061 
\ 1.006 
if 1.0059 



0.07 
~ 0.068 
^ 0.066 
g 0.064 
113 0.062 
0.06 

2046.2 
2046.1 
E, 2046 
& 2045.9 
2045.8 



fc 3098.5 
Jc 3098.4 
^ 3098.3 
"-' 3098.2 
3098.1 

2 4 6 8 10 

image progressive number 

Fig. 6. — Interchip transformation parameters as obtained from individual images. Data 
points from F225W are indicated with filled- circles, F275W with triangles, and F336W with 
crosses. The averages are indicated with solid lines, while the dashed-lines give the formal 
uncertainties. 

4. Average Absolute Scale relative to ACS/WFC 

The final step is to link, for each filter, the WFC3/UVIS chip [1] to an absolute plate 
scale in mas. To this purpose we adopt an average plate scale for our ACS/WFC master frame 
of 49.7248 mas ACS/WFC-pixel -1 (from van der Marel et al. 2007), and multiplied it by 
the -measured- scale factor between the WFC3/UVIS chip [1] and the master frame (which 
is expressed in ACS/WFC pixels). The results for the individual images and the averages 
for each filter, are shown in Figure while Table [3] gives the average values in mas pixel -1 . 
We believe that the differences in the relative values for the three filters are significant. The 
fact that the plate scales correlate with the wavelength suggests that refraction introduced 
by either the filters, or by the two fused-silica windows of the dewar, could have some role. 

Concerning their absolute values, instead, we have to consider that the velocity of HST 
around the Earth (±7 km s" 1 ) causes light aberration which induces plate-scale variations 
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up to 5 parts in 100 000 (Cox & Gilliland 2002)|£|, and that our master frame (from Villanova 
et al. 2007) was not corrected for it. 

The ACS/WFC plate scale for the Anderson's (2006, 2007) GD solution -once corrected 
for the temporal variations of the linear terms- has proved to be stable at a level of accuracy 
better than these velocity aberration variations (van der Marel et al. 2007). However, since 
we are not attempting to correct for this effect on our adopted ACS/WFC master frame, 
we simply limit the accuracy of the here derived WFC3/UVIS plate-scale absolute values to 
these accuracies, i.e. 12 parts in 100 000. 

Table 3: Average absolute plate scale of WFC3/UVIS in mas pixel -1 . Accuracy is 12 parts 
in 100 000. 



[mas pixel" 1 ] F225W F275W F336W 



39.760 39.764 39.770 



5 If we sum to this the Earth velocity around the Sun the plate-scale variations can reach up to 12 parts 
in 100 000 (Cox & Gilliland 2002). 
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Fig. 7. — Absolute scales relative to the one adopted for our ACS/WFC master frame 
expressed in mas. Each panel shows individual images, for different filters. As in the previous 
figure, the average values are indicated with solid lines, while the dashed-lines give the formal 
uncertainties. For reference, we show with error-bars the on-orbit maximum systematic 
errors introduced by the velocity aberration on the plate scale, i.e. for a velocity of ±7 km 
s -1 . However, since our master frame has not been corrected for scale variation induced by 
velocity aberration, these values can not be considered more accurate -in absolute sense- 
than 12 parts in 100 000. 
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5. Conclusions 

By using a limited (but best available) number of exposures with large dithers, and an 
existing ACS /WFC astrometric flat field, we have found a set of third-order correction coef- 
ficients to represent the geometric distortion of WFC3/UVIS in three broad-band ultraviolet 
filters. The solution was derived independently for each of its two CCDs. 

The use of these corrections removes the distortion over the entire area of each chip 
to an average accuracy of ~0.025 pixel (i.e. ~1 mas), the largest systematics being located 
in the ~200 pixels closest to the boundaries of the detectors (and never exceeding 0.06 
pixels) . We advise the use of the inner parts of the detectors for high-precision astrometry. 
The limitation that has prevented us from removing the distortion at an even higher level 
of accuracy is the lack of enough observations collected at different roll-angles and dithers 
which could enable us to perform an auto-calibration. 

Nevertheless, the comparison of the mid-2002 ACS/WFC positions with the new WFC3 
observations corrected with our astrometric solutions are good enough to clearly show the 
internal motions of u Centauri. These proved to be in perfect agreement with the most 
recent determinations. 

We also derived the average absolute scale of the detector with an accuracy limited 
by the uncertainties in the plate-scale variations induced by the velocity aberration of the 
telescope motion in the Earth-Sun system. 

For the future, more data with a longer time-baseline are needed to better characterize 
the GD stability of HST WFC3/UVIS detectors in the medium and long term. 
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